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The flow of a two-layer stratified fiuid over an abrupt topographic obstacle, simulating relevant 
situations in oceanographic problems, is investigated numerically and experimentally in a simplified 
two-dimensional situation. Experimental results and numerical simulations are presented at low 
Froude numbers in a two-layer stratified flow and for two abrupt obstacles, semi-cylindrical and 
prismatic. We flnd four different regimes of the flow immediately past the obstacles: sub-critical 
(I), internal hydraulic jump (II), Kelvin-Helmholtz at the interface (III) and shedding of billows 
(IV). The critical condition for delimiting the experiments is obtained using the hydraulic theory. 
Moreover, the dependence of the critical Froude number on the geometry of the obstacle are in- 
vestigated. The transition from regime III to regime IV is explained with a theoretical stability 
analysis. The results from the stability analysis are confirmed with the DPIV measurements. In 
regime (IV), when the velocity upstream is large enough, we find that Kelvin-Helmhotz instability 
of the jet produces shedding of billows. Important differences with flows like Von Karman's street 
are explained. Remarkable agreement between the experimental results and numerical simulations 
are obtained. 

PACS numbers: 47.55.Hd, 47.20.Ft 



I. INTRODUCTION 

The interaction between topography and flows with 
density interfaces is of interest both for fundamental rea- 
sons and for its relevance to practical problems. A most 
challenging problem is to describe and quantify the char- 
acteristic features occurring at the interface of a strat- 
ified flow passing over an obstacle. The interaction be- 
tween stable stratified flows and obstacles is a widespread 
phenomenon in nature. In the atmosphere, for example, 
the flow around buildings or mountains is particularly 
important because such conditions are often associated 
with high levels of atmospheric pollution due to low wind 
speeds and suppressed vertical mixing [H, Q • In physical 
oceanography, the interaction of marine currents with to- 
pographic features, such as ocean banks and coastlines, 
results in a complex system of circulation 0, S 0, S B 0] ■ 
In this case, observational, anal ytic al, numerical, and 
previous laboratorj^tudiesJEi, M El, [H, [H, [3, [H, 
H 113, [H, m, & M m [Mil] suggest that the com- 
bination of streamlines splitting, current intensification, 
and breaking of internal lee waves, play a significant role 
as a mixing source in the ocean [25|. Also, seamounts 
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have a decisive impact on enhancing the biological pro- 
ductivity and acting over the ecological processes that de- 
termine the structure of local ocean life [1^, . The rich 
diversity of geo-biophysical scenarios entirely justify the 
efforts for understanding the physics of the local insta- 
bility at the density interface, along with the conditions 
for which the generated lee waves break down. 

Among the topographic effects on fiows a particular 
problem that has received considerable amount of atten- 
tion is the oceanic observations in the Knight Inlet, a 
fjord in central British Columbia, Canada. The Knight 
Inlet experimentprovided field observations of tidal fiow 
over a large sill [J] ■ This work led to a series of papers re- 
porting different issues and possible explanation of them 
28, 29, 30, 31]. This fjord has a strong tidal flow which 
generates internal waves propagating along a pycnocline 
on both sides of a sill there. Farmer and Armi [3] have 
observed that a very large lee wave is formed behind the 
sill, while a train of strong depression pulses are gener- 
ated on the upstream side. 

Efforts has been dedicated to the upstream generation 
of solitary waves [2l|, [s^, |33] and also on the trapped 
wedge of mixed fluid behind the sill. Upstream influ- 
ence as a consequence of variable forcing, has been widely 
discussed in theoretical analysis, laboratory experiments 
and numerical simulations of stratified flow 28, 33]. Flow 
separation by a topographic obstacle and the final stage 
of vortex shedding has also received attention [2^. [291134]. 
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Model simulations of the lee wave formation by Lamb [3l| 
exhibit a stable flow over the sill until the lee wave over- 
turns. All those work agree that the larger-scale response 
can be sensitively dependent on small-scale instabilities. 
In the present work we focus on small-scale instabili- 
ties, addressing an analysis on the different regimes that 
appear downstream flow, near the obstacle. The impor- 
tance of the geometry of the obstacle in the transition 
from subcritical to supercritical will be demonstrated. 
Finally, we prove the Kelvin-Helmholtz (KH) instability 
at the two layers interface, is the triggering mechanism 
for the vortex shredding. 

Under suitable conditions, mainstream ceases to flow 
approximately parallel to the obstacle beyond a certain 
point (as the flow in a di verg ent channel) causing the phe- 
nomenon of separation [3ll Issj. In these cases, the flow 
next to the wall usually reverses its direction in some 
regions. Analytic theories for lee waves and hydraulic 
jumps usually assume that flow separation does not oc- 
cur However, when the topography is abrupt, that 
is, the vertical component of the velocity can not be ne- 
glected, flow separation is expected. Thus, in the present 
paper we compare experimental and numerical results 
with those obtained from hydraulic theories in the case 
of abrupt obstacles. 

The interaction between stratified flows and topogra- 
phy has been also investigated by means of laboratory 
experiments. In Ref. [l6| the lee-wave breaking pro- 
cess which occurs at low Froude numbers has been stud- 
ied in the case of a strong linear stratification and two- 
dimensional smooth (Gaussian-shaped) obstacles. In a 
slightly different setup, the structure of the far- wake vor- 
tices generated by a moving sphere in a linearly stratified 
medium was investigated in Ref. [12] ■ In this work the 
authors suggest that these vortices exhibit universal fea- 
tures similar to large-scale vortices found in the ocean. 
The case in which both stratification and rotation are im- 
portant has been study in Refs. where the role of 
laboratory experiments performed in close relation with 
numerical simulations is emphasized. 

In our previous work [36| we studied numerically and 
experimentally the different instabilities developed in a 
two-layer stratified flow over a pronounced obstacle. The 
existence of a concentrated jet in the lee side of the ob- 
stacle was demonstrated and we showed that the Kelvin- 
Helmholtz (KH) instability at the interface constitutes a 
more efficient source of mixing between the two-layers 
than the internal hydraulic jump. Here, we extend 
our work to two-layer fiows over prismatic and semi- 
cylindrical abrupt obstacles focusing on the role of the 
geometry, the appearance of secondary instabilities after 
KH, and the shedding of billows. 

The experiments were performed in a way that the flow 
could be assumed bidimensional. We use two layers of 
different density, with constant density inside each layer. 
We look for different regimes of the flow past the obsta- 
cles. The experimental results, obtained via direct visu- 
alization and also with Digital Particle Image Velocime- 



try (DPIV), showed to be in qualitative good agreement 
with numerical simulations and theoretical results. Four 
different regimes were found, depending on the global 
Froude number Fq and aspect ratio Hm- A complete 
diagram in the parameter space was obtained. The tran- 
sition from the regime HI to IV was explained with a 
stability analysis of the jet past the obstacle. Those an- 
alytical results were validated by velocity measurements 
with DPIV technique. 

This paper is organized as follows. In Section [H] we 
briefly review the treatment of two-layer flows using hy- 
draulic theory. In Section|TTT]we present the experimental 
setup. The numerical simulations are given in Sec. IIVI 
The results obtained from the experiments and simula- 
tions and their comparison are given in Sec. |Vl In Sec- 
tion |Vl] we present the results of linear stability analysis, 
which are relevant for the flow past the obstacle. Finally, 
a summary and the conclusions are given in section IVIII 

II. HYDRAULIC THEORY OF A TWO-LAYER 
FLOW PROBLEM 

Let us consider a linear theory for a flow consisting 
of two layers of different densities pi and p2 {p2 > Pi) 
over a fixed obstacle. Through this paper the subindexes 
1 and 2 correspond respectively to the upper and lower 
layer. We define a Cartesian reference frame with co- 
ordinates {x,y,z), where the fiow is in the x direction 
and z is directed vertically upwards. As we mentioned in 
the introduction, we shall focus on the situations where 
variations in the y direction can be neglected and the 
problem can be considered bi-dimensional. The depths 
of the layers are functions of the thickness of them 
depends on the position where it is measured. Thus, the 
depth of the upper and lower layers are di{x) and d2{x), 
respectively, and the height of the obstacle is h{x) as 
sketched in figure [T] Far upstream the depths are named 
di{x = 0) = diQ and d2{x = 0) = c?2o- 

The mean velocities of the fluid in each layer are Ui 
and U2. Let us assume that the fluid velocity is uniform 
far upstream, with uio — U20 — U. In addition, we 
assume: a) the pressure is hydrostatic, b) the Boussinesq 
approximation, which implies that e — (p2 ^ Pi)/ P2 ^ 1, 
is valid, and c) the top boundary of the upper layer is a 
free surface at constant pressure Ps, taken to be = 0. 
The flow over the topography is then characterized by the 
densities pi, p2^ the depth of the layers di, ^2, the mean 
velocity in each layer ui , U2 and the height of the obstacle 
h. With these assumptions, the Bernoulli functions for 
each layer may be written as: 

El = pig{d2 + di + h) + ^Piuj 

E2 = pigdi + P2g{d2 +h) + ^P2ul (1) 

Following Lawrence's model [13, [11], from the condi- 
tions dEi /dx = and imposing mass conservation in each 
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layer the following relation is obtained 

(1 - F^) dP ^ dh 
eF^F^ dx dx ^ ' 

where — F^ + F| — eF^^F| is the composite inter- 
nal Froude number while Ff = uf/{g'di) {i = 1,2) are 
the Froude numbers for each layer with g' — ^)g, 
and D = di + d2 + h. It has been shown [s^ that F is 
the adequate composite Froude number for characteriz- 
ing a two-layer flow; i.e. if F > 1 the flow is internally 
supercritical (the internal small waves cannot propagate 
upstream against the background flow), and if F < 1 
the flow is subcritical (the disturbances may propagate 
in both directions). When F = 1, the flow is termed 
critical and this location is usually called a control point. 
The transition from subcritical to supercritical flow is of 
special interest in our experiment. When the flow is su- 
percritical, internal hydraulic jump may take place which 
is an important source of turbulence and mixing. From 
Eq. (12), it follows that the critical condition F = 1, may 
occur \i dh/dx — 0. For that value of F the flow under 
goes a transition from subcritical to supercritical when 
the surface has horizontal tangent, i.e. at the crest of the 
obstacle. If, in addition, the condition e ^ 1 is imposed 
(which is satisfied in our experiments), the composite 
Froude number may be expressed as F^ = F^ + F| . On 
the other hand, in the present study we consider flows 
where Ff < 1. Thus, from Eq. ^ and the Boussinesq 
approximation it follows that 

dD ^ eF^Fj dh_ 

dx 1-F^ dx'^ ' ^ ' 

therefore we consider D = dio + c?20 a- constant. Hence 
the free surface will be taken as horizontal. 

In order to obtain the critical values of the flow pa- 
rameters for the subcritical-supercritical transition, we 
impose that E2 — Ei is constant and using the Boussi- 
nesq approximation, then 

(4) 

where r = d2/D, tq — d2o/F>, H = h/d2a and 

= (5) 

V 5 "10 9 0-20 

is the Froude number F calculated in the upstream flow 
far from the obstacle. We shall call Fq as the global 
Froude number. On the other hand, the critical condition 
F^ — 1 may be expressed as 

Thus, Eqs. (HI) and imposed at the crest of the 



obstacle, may be used to determine, Fqc, the critical value 
of Fq for the occurrence of supercritical flow. In our 
experiments we fixed the aspect ratio tq = 0.6, and we 
considered different values of Hm = hm/d2Q, hm being 
the height of the obstacle. The critical values, Fqc, for 
the occurrence of critical flows as a function of Hm are 
obtained solving equations ([U and (O. The results are 
shown in Sec. |Vl Fig. [71 where are compared with the 
experimental results. 

Theoretical solutions for the flow when Fq is larger 
than Fqc, i.e. when the flow is beyond the critical condi- 
tion, have been obtained in Refs.jll, H^. These solutions 
predict that there is range of values of Fq for which there 
is a wave that moves backwards to the flow changing the 
conditions upstream. These approaches allow the calcu- 
lation of the velocity and amplitude of this wave. This 
kind of wave has been also obtained using the KdV equa- 
tion [39], and numerical simulations, with step like strat- 
ification 32] and linearly varying stratification j21|]. This 
effect may difficult the experiments inside the container 
in some cases. However, we show that quantitative good 
results may be obtained within a broad range of param- 
eter values. 



III. EXPERIMENTAL SETUP 

Our experiments were performed in a water tank of 
size equal to 2.0 x 0.29 x 0.137 m'^. We used a closed 
channel where we towed the obstacle at different veloc- 
ities with a calibrated motor (Fig. [21). The velocities of 
the fluid were calculated in a reference frame fixed to the 
obstacle. In our previous work [36| . we showed that this 
configuration is equivalent to the flow over a fixed obsta- 
cle. Moreover, we focused our observations only on the 
central region of the tank where the structures are per- 
sistent and the far boundaries effects can be neglected. 
In this work, we always show images and diagrams of 
leftward moving obstacles. We used two different obsta- 
cle shapes, prismatic and semi-cylindrical, both of them 
having a height of hm = 0.125 m, width Wq = 0.13 m 
and length of L — 0.25 m, see Fig. [31 The obstacles were 
scaled in such a way that the confinement aspect ratio 
verifies Wq/W < 1 {W = 0.137 m, width of the water 
tank) and the lateral flow around the obstable can be 
neglected in order to reduce the problem to a quasi-two 
dimensional situation. As mentioned above, we want to 
model a density profile with an abrupt gradient at the in- 
terface. In order to get this step- like stratification we first 
filled the tank with a layer of density p2 = 1002 kg/m!^ 
using NaCl solution. To fill the upper layer, pure water 
with density pi = 1000 kg/m!^ was carefully poured over 
a sponge floating on the free surface. The time scale of 
the molecular diffusion between the layers is much longer 
than the typical experimental times. Despite that, due 
to the mixing produced by the moving obstacle, after a 
few measurements, the tank had to be emptied and new 
fluid layers poured again. 
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The flow was examined via two standard techniques: 
dye visuahzation and digital particle image velocimetry 
(DPIV). In the first case, the upper layer was dyed with a 
KMn04 solution in order to obtain a good visual contrast 
between both layers. A powerful source of fluorescent 
light from behind was used to obtain a uniform illumina- 
tion. We obtained global qualitative picture of the flow 
with this technique. In the second case, DPIV allows us 
to obtain quantitative values of velocity field based on 
the cross-correlation of two consecutive images recorded 
by a digital camera. Neutrally buoyant polyamide parti- 
cles of 50 X 10~^ m diameter were seeded in the bottom 
layer. A green laser sheet of 100 mW is used to illuminate 
a cross section plane of the flow which is recorded by a 
digital camera PIXELINK PL-A741. In order to control 
the interface and avoid attenuation of the LASER light as 
much as possible, when using this technique, only a thin 
layer of water at the bottom of the upper layer is dyed. 
For both obstacle shapes, five different sets of heights 
were chosen, always keeping constant the characteristic 
ratio To = 0.6. The experiment was repeated for each 
set of heights with a wide range of velocities analyzing 
the different behaviors, using the two mentioned visual- 
ization techniques. 



IV. NUMERICAL METHOD 

The numerical simulations considered here were ob- 
tained with the in-house flow solver caffaSd.MB [i^] de- 
veloped jointly by Universitat Rovira i Virgili (Tarrag- 
ona, Spain) and Universidad de la Republica (Montev- 
ideo, Uruguay). It is an original Fortran95 implementa- 
tion of a fully implicit finite volume method for solving 
the 3D incompressible Navier-Stokes equations in com- 
plex geometry using block structured grids. This three- 
dimensional solver, based on a previous two-dimensional 
solver f4l!| , is described and validated in [i^, ■ 

The unsteady incompressible Navier-Stokes equations 
with Boussinesq approximation for buoyancy terms were 
considered. Since the Reynolds number was below 
Re=8.10^ for all cases, no turbulence model was required, 
thus transient solutions were computed directly. The 
time step was set to 2.0 x 10~^ s for all cases. This time 
scale is about (hm/U) /lO'^ for the highest velocity case. 
Simulations were run starting from null velocity fields 
through lO"* time steps, or about 200 s of flow time. 

In the simulations, the obstacle remains fixed against 
a steady two-layer current of fluid. Thus a uniform ve- 
locity profile was specified at the upstream boundary lo- 
cated at a distance of 8hm upstream from the obstacle, 
and a null gradient outlet was used at the downstream 
boundary, located 15hm downstream. As the top surface 
is not disturbed by the flow, it was modeled as an hor- 
izontal slip boundary at fixed height. All other bound- 
aries correspond to wall surfaces and non-slip condition 
was directly applied to them, including the vertical walls 
of the channel. Thus, the simulation accounts fully for 



three-dimensional effects, although the flow reveals itself 
as essentially two dimensional owing to the geometry of 
the obstacles and the relatively low Reynolds numbers. 

For both obstacles the grid was made up of three 
blocks, although the topology was different in each case. 
For the prismatic obstacle three straight blocks were as- 
sembled, two at each side of the obstacle and the third 
extending on top of them along the domain. On the 
other hand, for the cylindrical obstacle one C-grid block 
was used around the obstacle together with two other 
straight blocks, upstream and downstream of the obsta- 
cle. 

Grid resolution was set essentially uniform through the 
domain at /im/25, being enough to resolve flow details at 
these rather low Reynolds numbers. Due to the layout 
of the grid in the cylindrical obstacle case the spatial 
resolution normal to the wall was slightly higher near 
the obstacle, reaching about h„i/35. 



V. RESULTS 

In this section we are going to discuss the experimental 
and numerical results, which exhibit a great qualitative 
agreement between them. It is also worth noting that, 
in spite of quantitative differences, the results for both 
topographic shapes show a qualitative similarity. Let us 
start discussing the different regimes observed in the ex- 
periments. Using the dye technique we visualized the dif- 
ferent regimes as a function of the obstacle velocity and 
the aspect ratio. In Fig. 2] (Fig. [5]), wc show experimental 
and numerical results for the prismatic (semi-cylindrical) 
obstacle. The aspect ratio and the height of the bottom 
layer are the same for both obstacles, i.e. rg — 0.6 and 
d2o = 15 cm. 

In these figures we distinguish four regimes, all of them 
present in both obstacles. However, both, experimental 
and numerical results reveals one very important point: 
the regimes take place at different critical Froude number 
for different geometries. For the prismatic obstacle, the 
velocities vary between 0.12 cm/s {Fq = 0.035) and 0.64 
cm/s {Fq = 0.187). While for the semi-cylindrical obsta- 
cle, the velocities vary between 0.13 cm/s (Fq = 0.039) 
and 0.43 cm/s {Fq = 0.125). The numerical simulation 
for both obstacles were performed at the same global 
Froude numbers as the experimental results. 

All the different behaviors of the downstream flow cor- 
respond to the case in which the flow upstream is sub- 
critical, i.e. Fq < 1. The thickness of the layers varies 
from one point to another immediately after the obsta- 
cle, as can be seen in Figs. mandO and consequently the 
Froude number F too. We look for the first critical value 
of the local Froude number {F = 1), over the obstacle 
in a place where ^2 is minimal and where the velocity is 
higher. 

The regime (I) corresponds to the situations in which 
the flow is subcritical, F < 1, along all the flow over the 
obstacle. At low velocity, we observe that a jet through 
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the bottom layer is flowing near the obstacle. As the 
velocity is increased, the jet begins to separate from the 
obstacle and rise towards the horizontal. 

A transition from regime (I) to (II) occurs when Froude 
number reaches a control point {F — 1) somewhere over 
the obstacle. In regime (II) the flow is subcritical up- 
stream, and supercritical past the obstacle. As a conse- 
quence, an internal hydraulic transition is developed at 
the lee side. The interface between the layers is smoothly 
disturbed both in regime (I) and (II), however, over the 
obstacle the interface is stable enough not to break but 
to induce a jet in the lower layer. When the velocity gra- 
dient between the jet and the surrounding fluid is strong 
enough a Kelvin-Helmholtz instability appears inside the 
lower layer at the lee side of the obstacle. The wavelength 
of KH depends on two densities and the velocity profile 
(see Sec.lVip. The regimes (II) and on, are all supercrit- 
ical past the obstacle, F > 1 beyond the critical point. 

Further increasing the Froude number Fq (still below 
1), we reach regime (III) where it is clearly visible a lee 
wave that perturbs the interface separating the two lay- 
ers, with a quasi-sinusoidal profile. Downstream, a sec- 
ondary instability develops and a kind of mixing is ob- 
served. Finally, the regime (IV) is characterized by the 
shedding of vortical portions of lighter fluid separated 
from the upper layer. As a consequence, intense mixing 
between the two layers takes place. In this regime, the 
interface between both layers is strongly disturbed. The 
jet drags fluid from the upper layer and billows formation 
appear. It is worth mentioning that the frequency of the 
shedding is almost constant with respect of variation in 
the velocity U. 

The numerical results are in good agreement with the 
experimental results, as can be visualized in Figs. H] and 
[SJ The regimes observed in the experiments were also 
obtained in the numerical simulations for the two obsta- 
cles that we used. As Fq gets larger the behavior of the 
numerical and the experimental results differs. For the 
case of the prismatic obstacle the vortex shedding starts 
before than in the simulations (Fig.|4tl). Whereas in the 
cylinder case the numerical experiments shows a vortex 
shedding before the experiments (Fig.[5]d). These differ- 
ences will be explained further on. The overall agreement 
between numerical and experimental method is good as 
can be seen in Fig. [S] which compares numerical and ex- 
perimental results corresponding to regime III and IV. 
The flow separation occurs before in the experiments for 
low Fq (Fig. [6]a-b). However for greater Fq the situation 
inverts (Fig. [S] c-d) . 

In Fig. ([7]) we show the global Froude number Fq as 
a function of the aspect ratio Hm for the transition from 
sub-critical to super-critical. The continuous line cor- 
responds to the theoretical global Froude number for a 
smooth obstacle as it was calculated solving equations 
(HI and ([6]). The experimental Fq is represented with 
open squares for the semi-cylindrical obstacle and filled 
squares for the prismatic one. The dashed lines are the 
linear fitting of the experimental points. This graph 



clearly reflects the fact that the smoother the obstacle, 
the higher the velocity that will be needed to reach the 
critical Froude number. Then, the critical values for the 
semi-cylindrical obstacle are aligned between the values 
of the prismatic obstacle and the smoother theoretical 
prediction. In this way, by means of the laboratory ex- 
periments we prove that he transition from one regime 
to another strongly depends on the geometry, the sub- 
critical flow may correspond to a supercritical flow for 
another geometry. 

In order to obtain the transition between the differ- 
ent regimes, we use the DPIV technique. For the same 
parameters of the previous situation, c?20 = 15 cm and 
prismatic obstacle, we observed that when the velocity 
is lower than U = 0.11 cm/s, we obtained a sub-critical 
flow. The transition between sub-critical to super-critical 
flow takes place at C/ = 0.12 cm/s {Fq = 0.035). For 
U = 0.38 cm/s (Fq = 0.111), KH instability at the inter- 
face appears. Finally, for velocities higher than U = 0.64 
cm/s {Fq = 0.19) shedding of billows takes place. Fig- 
ure [5] summarizes in a stability diagram the different 
regimes as a function of the global Froude number Fq and 
the aspect ratio Hm for tq = 0.6. The diagram shows the 
transition values between the regimes for the prismatic 
(full symbols) and semi-cylindrical (open symbols) ob- 
stacles. The dotted lines are linear approximations. We 
show for Hm = 0.83, in open symbols the transitions 
between the different regimes for the cylindrical obsta- 
cle. It is interesting to note that, although the prismatic 
obstacle reaches the super critical regime at a lower ve- 
locity than the semi-cylindrical obstacle, KH develops at 
higher velocities than in the case of the semi-cylindrical 
obstacle. 

The slope of the jet after the obstacle plays an impor- 
tant role in the transition between the different regimes. 
Indeed, it was observed that transition depends not only 
on Fq, tq and Hm but also in the shape of the obstacle, i. 
e. the slope of the jet after the obstacle. As a result, the 
lee waves past the prismatic and semi-cylindrical obsta- 
cles have similar characteristics at different Froude num- 
bers. We exemplified this in figure[9l We show KH insta- 
bility at the interface for the two obstacles. Despite the 
prismatic obstacle reaching the super critical regime at 
a lower velocity than the semi-cylindrical obstacle, the 
abrupt geometry of the bank makes difficult to rise up 
the jet and develop KH. Thus, for the prismatic obsta- 
cle, KH develops at higher velocities than in the case of 
the semi-cylindrical obstacle. 

We characterized the velocity profile via DPIV analy- 
sis. Figure [TO] (b) shows the visualization of the velocity 
profile over the semi-cylindrical obstacle obtained with 
DPIV technique at the places indicated in Fig. [TO] (a). We 
observe that the shape of the jet is nearly triangular and 
quite constant at different distances from the obstacle. 
These facts will be used to make an hydraulic stability 
analysis of the interface (see Sec. IVI|) . 

In regime HI as the flow velocity U is increased, the 
inclination of the jet, the amplitude of the interface wave 
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U (cm/s) a (cm) xi, (cm) Xb/a 

0.45 2A 5.4 2.6 

0.55 2^3 6.2 2.7 

0.65 2.8 7.4 2.6 



TABLE I: Thickness of the jet a and distance between the 
billows xi, for ro = 0.6 and Hm ~ 0.83. 



interface taking into account the effects of the upstream 
wave. On the other hand, for Fq — 0.11, d'lQ/diQ — 1.03, 
this estimation shows that the effects of this wave are 
negligible. Furthermore we note that the results in Fig.[7| 
and Fig. [8] for Fq < 0.11 are not distorted in any way by 
the upstream effect. In the region Fq > 0.11 altough the 
results could be slightly distorted, they remains qualita- 
tively correct as we can conclude from comparison with 
the numerical simulations. 



and the velocity of the jet increase. However, when the 
Froude number exceeds the critical value, the wave at 
the interface is no longer observed and the flow enters in 
regime IV. In this regime, there is shedding of vortical 
portions of the more lighter fluid which submerge into 
the lower layer. 

The shedding of billows in regimen (IV), as can be seen 
in figure [TTl provide evidence that the instabilities which 
lead to this regimen could be caused by a KH instability. 
In order to test this hypothesis, we correlate the thickness 
of the jet a and the distance between the billows Xb- If the 
shedding was caused by KH instability, the distance be- 
tween the billows must be approximately the wavelength 
of the most unstable mode. Table I shows the values of 
a, Xb and the ratio between them for the case rg = 0.6 
and Hra = 0.83, for the prismatic obstacle at different 
velocities. These results reveals a clear proportionality 
between a and Xb, with a slope of Xb/a = 2.7. Indeed, in 
the next section we consider a triangular jet, and we show 
that the dimensionless wavenumber of the most unstable 
mode is given by kb — 1.225. Then, the ratio between 
the corresponding wavelength and the thickness of the 
jet a {a = 2b) is given by A/a = 2.57, which is in very 
good agreement with the ratio (2.6 — 2.7), obtained ex- 
perimentally. 

As commented in Sec. [Tl] theoretical solutions based 
on hydraulic theory predicts that for Fq > Fqc there is 
a wave that propagates in opposite sense to the flow, 
modifying the upstream conditions. The role of this up- 
sream wave has been profusely studied theoretical, ob- 
servational and numerically [j, IH, [s^, |3lj . We address 
here, the question of how much this feature affects our 
results. 

In order to study the effect of this upstream wave, we 
performed numerical simulations, in a broad range of val- 
ues of the Fq, focusing on the dynamics far upstream the 
obstacle and looking for the transitory effects. The tem- 
poral evolution of the upstream wave, for Fq — 0.35, 
is shown in Fig. [121 We note that the profile of the up- 
stream wave is very smooth in comparison with the down- 
stream perturbation. Accumulation of wave effects in the 
closed container for Froude numbers clearly exceeding the 
critical condition produces that the thickness of the lay- 
ers are larger than those obtained in an unlimited flow. 
This effect is more important as larger the Fq is. For 
instance, from the theory [l| follows that for Fq — 0.181 
and Hm = 0.833, the increment of the thickness is signif- 
icant, c?']^o/^io = 1-10, where d'lQ is the thickness of the 



VI. STABILITY ANALYSIS OF THE JET PAST 
THE OBSTACLE 

A. Kelvin-Helmholtz Instability 

As described in section |Vl in the lee side of the ob- 
stacle the flow presents different regimes depending on 
the values of Fq, tq and Hm- However the type of flow 
downstream is not only determined by the flow param- 
eters but also by the slope of the jet past the obstacle. 
We will analyze the stability of this flow. 

In regime HI, the interface between the two layers is 
perturbed by a quasi-sinusoidal wave, which breaks down 
when its amplitude is sufficiently large, indicating that a 
secondary instability takes place. The DPIV measure- 
ments clearly show that the jet is located near the inter- 
face, we conjecture that this wave results from a primary 
Kelvin-Helmholtz instability. We test this hypothesis 
based in the experimental measures shown in Fig. [TU] 
(b) and we will describe the flow near the interface as a 
triangular jet. The flow is given by 

^ I Vi for z > 6 and z < —6, 

~ \V2 + {Vi-V2)\z\/b fOT-b<Z<b, 

where 2b is the thickness of the jet and the densities are 
P2 for z < b and pi for z > b. We perform standard 
linear stability analysis assuming a stream function of 
the perturbations i/j of the form j44| 

ip{x, z,t) — (f){z) e:iq3{ik{x — ct)) , (8) 

where exp(ik{x — ct)) = exjp{ik{x — c,.i)) exp at, with a — 
kci, and Ci,Cr are the imaginary and real parts of the 
phase velocity c. Thus, a represents the growing rate of 
the perturbations. If ct > 0, then the flow is unstable to 
these perturbations, and if ct < the flow is stable. Since 
the flow under consideration has constant density in each 
layer, 4'{z) is determined by Rayleigh's equation: 

r{z)~k'm^O, (9) 

whose general solution is V'i(^) — Aiexp{kz) + 
Si exp(— fcz), where the index i = 1,2,3,4 denotes the 
solution for the regions z>b,b>z>Q,0>z>~b and 
z < —b respectively. The perturbations must be bounded 
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at z ^ oo and z —>■ — cx), then the (f>i^s become 



hiz) 

h{z) 
h{z) 
h{z) 



Bi cxp(— 

A2 exp(/cz) + B2 exp(— 
^3 exp(A:z) + exp(— /cz) 
A/^ exp(A:2;). 



(10) 



We now impose the boundary conditions across the 
Umits separating the regions by requiring continuity of 
vertical velocity and pressure. Since the velocity profile 
is continuous, the first of these conditions is equivalent 
to require (j)i{zj) = (j)i+i{zj) [ij, where Zj denotes the 
position of the boundary that separates regions i and 
i + 1. The second of these conditions requires 

{Uiizj)-c)<j)'{zj),-U-{zj)(t)i = {Ui+i{zj)-c)<l)'{zj),+i~U-{zj 

The boundary condition yield a linear system of equa- 
tion on {Bi, A2, B2, A3, B3, A4). The dispersion relation, 
i. e. the relationship between k and c, is obtained banish- 
ing the trivial solution. The numerical results obtained 
by solving the dispersion relation are shown in figure 
(fT3)) . The normal mode stability analysis yields that 
the more unstable mode corresponds to the wavenum- 
ber k = 1.225/5. For b = 1.1 cm (see Fig. [10] b) , the 
corresponding wavelength is A = 10.3 cm which is in good 
agreement with the experimental data of A = 10.2 cm. 
We also performed the linear stability analysis of a jet 
that is located near the wall of the obstacle, with field 
velocity given by 

y-l^' > (11) 

\V2 + {Vi-V2)\z\/b for-fe<z<6, ^ ^ 

where the rigid boundary is located at z ^ —b, and the 
densities are p2 for z < b and pi for 2; > 5. In this case, we 
have three regions with corresponding stream functions 
(/)i(z), (/)2(z), </)3(2:). The condition at the rigid boundary 
is given by (ji^^—b) = and the other boundary 

conditions are the same as the ones imposed for the un- 
bounded jet. From this analysis we obtained that the jet 
is partially stabilized in the proximity of the wall of the 
obstacle, that is, the instabilities grow slowly compared 
with the unbounded case (Fig. (|13p). This result is in 
agreement with those obtained previously by Hazel [45l |. 
This stabilizing effect is caused by the combination of the 
density stratification and the presence of the wall. Then 
the wall has little effect on the stability of the jet, if the 
density is constant. This explains why the instabilities 
appear downstream, but not over the obstacle, where the 
jet is located on a rigid boundary. 

The stability analysis results and the measured wave- 
length of the lee waves in the regimes III and IV, 
strengthen the assumption that the instabilities observed 
past the obstacle in both cases are caused by Kelvin- 
Helmholtz instability. We can also explain the fact that 
the frequency of the billows is almost independent of the 
velocity U . From Table I we can see that the thickness of 



the jet is (for fixed tq and Hm) practically proportional 
to U . The frequency of the billows can be expressed as 
f — Xb/ub , where Ub is the velocity of the billows. Since 
Xb is nearly proportional to U, if we assume that Ub is 
proportional to U, then the frequency is weakly depen- 
dent on the velocity U. This is in contrast with vortices 
produced in flows like the Von Karman's street, where 
the frequency, in a wide velocity range, is almost pro- 
portional to the velocity. The difference resides in the 
fact that in the Von Karman's flow, the characteristic 
length is constant (the diameter of the cylinder), while 
here the characteristic length (the thickness of the jet) 
varies almost proportional with U, which causes that the 
wavelength of the most dangerous mode increases with 
U. This effect compensates the increasing of Ub with U, 
■j^avlng the frequency almost constant. 



VII. SUMMARY AND CONCLUSIONS 

In this work, we studied two-layer stratifled flow over 
abrupt obstacles, more specifically prismatic and semi- 
cylindrical. It is remarkable that in all cases, the critical 
value of the global Froude number, Fq, is less than that 
predicted by the hydraulic theory. In those cases, the ver- 
tical components of the velocity over the obstacles cannot 
be neglected. As a consequence of its smoother geome- 
try, the results for the semi-cylindrical obstacle are closer 
to the hydraulic theory than those of the prismatic one. 
This result indicates that an abrupt obstacle reaches the 
control point for lower velocities than for a smooth ob- 
stacle. In other words, for a given Fq, a flow subcritical 
for a geometry could be supercritical for another geome- 
try. Despite the fact of being approximate, the hydraulic 
theory allow us to estimate the thickness and averaged 
velocity of the strong jet that forms past the obstacle. 
However the presence of flow separation and dissipation 
at the lee side of the obstacle prevent the use of hydraulic 
theory for making further predictions. 

For the downstream flow, four different regimes were 
identifled and represented in the parameter space Fq — 
Hjn- We showed typical images of them in Fig. U) In 
all the regimes observed the flow upstream remains sub- 
critical. In regime I, the flow is subcritical everywhere, 
i.e. upstream and downstream. In regime II and all the 
subsequent regimes the flow past the obstacle is super- 
critical. The regime II is characterized by the transi- 
tion from subcritical to supercritical by the local Froude 
number experimentally measured. Also in regime II, a 
shear instability develop mainly in the lower layer. The 
interface that separates the layers is disturbed with a 
quasi-sinusoidal wave in regime III. This disturbance is 
due to Kelvin-Helmholtz instability. Furthermore, this 
quasi-sinusoidal wave disturbance breaks down through a 
secondary instability for sufficiently large amplitude and 
departure from two-dimensionality takes place. Regime 
IV is characterized by the shedding of billows. This pro- 
cess increases the rate of mixing between the two layers. 
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A stability analysis on the jet immediately after the 
obstacle allowed us to verify that the shedding was also 
produced by Kelvin-Helmholtz instability. The charac- 
teristic length of the most dangerous mode of the strat- 
ified jet is in agreement with the distance between the 
billows. We have also shown that the stabilizing effect re- 
sulting from the combination of the proximity to an rigid 
boundary and stratification turns out in not observing 
any instabilities over the obstacle. 

We studied the wave that propagates upstream when 
Fq > Fqc- In the simulations, we reproduced this wave 
with excellent agreement with the hydraulic theory. The 
enhancement of the deeper layer thickness as well as the 
velocity of the fluid inside it are in very good agreement 
with the theory [l| . These results indicate that the effect 
of this wave must be considered in flows exceeding the 
critical conditions. For Fq of the order 0.1 or less, the 
effects of the wave may be negligible. Furthermore, the 
results in Fig. [7] and Fig. [5] for Fq < 0.1 are not distorted 
in any way by the upstream effect. For this region of the 
space diagram, and high enough Hm, there are still the 
four different regimes, as it was described before and also 
in this region vortex shedding takes place. In the region 
Fq > 0.10 the results remains qualitatively correct as we 



can conclude from the numerical results comparison. 

The experimental and simulation results showed that 
the inclination of the jet plays a central role. Immedi- 
ately after the separation, as velocity U increases the in- 
clination of the jet decreases. However, when U is larger 
than a critical value, the inclination of the jet increases as 
the velocity U increases. When the flow enters in regime 
IV, the jet slope is pronounced. The different regimes 
reported in this work appeared for both obstacles, but, 
since separation is a phenomenon strongly dependent on 
the curvature of the surface, the critical values for passing 
from one regime to another depend on the geometry of 
the obstacle. 
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FIG. 3: Side (top) and plan view (bottom) of the prismatic 
and semi-cylindrical obstacles. Where dio, d20 are the depth 
of the layers, hm the maximum height and L the length of 
the obstacle. W is the width of the water tank and Wo the 
width of the obstacle. 



FIG. 1: Sketch showing variables used in the description of 
the two-layer flow. U is the upstream uniform fluid velocity, 
di, d2 are the depth of the layers and h is the height of the 
obstacle at a given location with hm the maximum height of 
it. 




FIG. 2: Schematic diagram of the experimental conflguration 
and the two visualization techniques used, a) The upper layer 
is dyed with KMn04 to obtain good contrast for visualization 
when the tank is lighted from behind, b) The DPIV is carried 
out lightening from above with a LASER. In this case, in order 
to distinguish the interface, only a thin portion of the upper 
layer is dyed. 
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Experimental Numerical 
Fo = 0.035 (a) 




Fo = 0.047 (b) 




Fo = 0.111 (c) 




Fo = 0.187 (d) 




25 cm 



FIG. 4: Experimental snapshots (left) and numerical simula- 
tions of a flow (right) past an prismatic obstacle corresponding 
to dio = 0.10 m (upper layer), d20 = 0.15 m (bottom layer) 
corresponding to the different regimes, a) Subcritical regime 
(I), b) Internal hydraulic transition (II). c) KH instability at 
the interface (III), d) Billow formation (IV). In all the im- 
ages the aspect ratio is ro — 0.6 and values of the the Froude 
numbers Fo are indicated. 
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Experimental Numerical 
Fo = 0.039 (a) 




Fo = 0.053 (b) 




Fo = 0.072 (c) 




Fo = 0.125 (d) 




25 cm 



FIG. 5: Experimental snapshots (left) and numerical simula- 
tions of a flow (right) past an semi-cylindrical obstacle corre- 
sponding to dio — 0.10 m (upper layer), d20 = 0.15 m (bottom 
layer) corresponding to the different regimes, a) Subcritical 
regime (I), b) Internal hydraulic transition (II). c) KH insta- 
bility at the interface (III), d) Billow formation (IV). Same 
parameters as in Fig|4] 
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Experimental Numerical 
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FIG. 6: Comparison between experimental (left) and numer- 
ical (right) results for semi-cylindrical obstacle. Top row 
d2o = 15 cm and U = 0.38 cm/s and Fo = 0.111. Bottom row 
d2o = 15 cm and U = 0.64 cm/s and Fo = 0.187. 
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FIG. 7: Critical Global Froude number as a function of the 
aspect ratio Hm = hm/d2o, for a fixed relation between the 
height of the layers ro = 0.6. The continuous line corresponds 
to the hydraulic model and the symbols to the experimental 
results of the two obstacles considered: semi-cylindrical (open 
squares) and prismatic (full squares). 



FIG. 8: Stability diagram of the different regimes for the pris- 
matic and the semi-cylindrical obstacle as a function of the 
global Froude number Fo and the aspect ratio Hm = hm/d2o- 
The relation between the height of the two layers is fixed to 
ro = 0.6. The full symbols correspond to the experimen- 
tally obtained transition between different regimes for the 
prismatic obstacle. Dotted lines are linear approximations 
of those experimental points. The open symbols correspond 
to the transition points for the semi-cylindrical obstacle for 
Hm ~ 0.83. Squares (full and open), correspond to the transi- 
tion between sub-critical to internal hydraulic jump regimen. 
Circles (full and open), correspond to the onset of Kelvin- 
Helmholtz regimen. Triangle symbol (full and open), corre- 
spond to the shedding of billows regimen. 




FIG. 9: Flow passing over the two obstacles of different ge- 
ometry. These snapshots correspond to the Kelvin-Helmholtz 
instability between both layers at U = 0.35 cm/s for the pris- 
matic obstacle and at U = 0.22 cm/s for the semi-cylindrical 
obstacle. 
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FIG. 10: a)DPIV image for the semi-cylindrical obstacle for 
KH instability at Fq = 0.122. The vertical lines correspond to 
the places where the velocity profiles were taken, b) Velocity 
profiles for the three positions marked in the image displayed 
in a) 
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FIG. 11: Shedding of billows for = h/d2 = 0.76, ro = 0.6 
and F — 0.26 for the prismatic obstacle. The frequency of 
shedding is 0.5 s^^. 
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FIG. 12: Transient evolution of the interface showing the up- 
stream influence of the wave corresponding to Fo = 0.35. The 
simulations start ai t = with a null- velocity condition. 



FIG. 13: Non-dimensional growing rate of the perturbations 
for the unbounded stratified jet (continuous line) and the 
stratified jet near a wall (dashed line), defined by Eqs.(l7| 
and (tni), where k* = kb and a* = ab/{V2 - Vi). 



